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ABSTRACT 

We present a photometric investigation of the variation in galaxy colour with envi- 
ronment in 11 X-ray luminous clusters at 0.07 < z < 0.16 taken from the Las Cam- 
panas/ AAT Rich Cluster Survey. We study the properties of the galaxy populations 
in individual clusters and take advantage of the homogeneity of the sample to combine 
the clusters together to investigate weaker trends in the composite sample. We find 
that modal colours of galaxies lying on the colour-magnitude relation in the clusters 
become bluer by d{B — R) / drp = — 0.022±0.004 from the cluster core out to a projected 
radius of Vp = 6Mpc; further out in radius than any previous study. We also examine 
the variation in modal galaxy colour with local galaxy density, S, for galaxies lying 
close to the colour-magnitude relation and find that the median colour shifts bluewards 
by d{B — i?)/(ilog]^o(^) = —0.076 ± 0.009 with decreasing local density across three 
orders of magnitude. We show that the position of the red envelope of galaxies in the 
colour-magnitude relation does not vary as a function of projected radius or density 
within the clusters, suggesting that the change in the modal colour results from an 
increasing fraction of bluer galaxies within the colour-magnitude relation, rather than 
a change in the colours of the whole population. We show that this shift in the colour- 
magnitude relations with projected radius and local-density is greater than expected 
from the changing morphological mix based on the local morphology-density relation. 
We therefore conclude that we are seeing a real change in the properties of galaxies on 
the colour-magnitude relation in the outskirts of clusters. The simplest interpretation 
of this result (and similar constraints in local clusters) is that an increasing fraction of 
galaxies in the lower density regions at large radii within clusters exhibit signatures of 
star formation in the recent past, signatures which are not seen in the evolved galaxies 
in the highest density regions. 

Key words: surveys, catalogues, galaxies: photometry, galaxies: clusters: general, 
galaxies: clusters: individual: A 22, A 550, A 1079, A 1084, A 1285, A 1437, A 1650, 
A 1651, A 1664, A 2055, A 3888. 



1 INTRODUCTION 

Visvanathan and Sandage (1977) first noted that the more 
luminous early-type galaxies (ellipticals and SOs) within 
Coma and eight other local clusters exhibit systematically 
redder integrated colours than the late-types. These galax- 
ies also demonstrate a tight correlation between their colours 
and magnitudes: the colour-magnitude relationship (CMR; 



e.g. Bower, Lucey and Ellis 1992). Kodama & Arimoto 
(1997) show that the slope of the colour- magnitude relation 
is primarily due to the mean stellar metallicity (as suggested 
by Dressier 1984) , while the scatter around the relation may 
be due to age effects (Kodama et al. 1999). 

As one moves from the cores of clusters into the sur- 
rounding field the morphological mix in the galaxy popu- 
lation shifts from almost completely passive, early-types to 
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one dominated by star-forming spiral galaxies. This is the 
origin of the Morphology-Density relation (T-E; Dressier 
1980). The environmental "transition" region between clus- 
ters and the field is an important region to study if we are 
to loarn what impact environment has on galaxy proper- 
ties such as star formation rate, luminosity and morphology. 
Recent numerical simulations of clusters of galaxies show 
that the history of galaxy accretion is retained in the clus- 
ter's radial profile. The first galaxies to be accreted still re- 
side nearer the centre while more recent additions arc lo- 
cated preferentially on the outskirts of clusters (Diaferio et 
al. 2001). This means that if the star formation in galax- 
ies is suppressed when they enter the cluster environment 
(Balogh, Navarro & Morris 2000), then the variation in the 
typical accretion age of cluster galaxies with radius may be 
visible in the ages of the stellar populations of the passive 
cluster members which inhabit the colour-magnitude rela- 
tion. Studies of galaxies in the outskirts of clusters and the 
variation in their properties with rcdsliift can therefore in- 
form us about whether environment has altered the charac- 
teristics of cluster galaxies as the clusters have been assem- 
bled over the Hubble time. 

There are two main observational techniques to trace 
the variation in the stellar populations of galaxies between 
the field and cluster and its evolution with redshift: photom- 
etry and spectroscopy. 

Wide-field photometric studies of clusters have un- 
covered evidence of changes in the colours of the galaxy 
populations with radius: Butcher & Oemler (1978, 1984) 
first identified a photometric radial gradient; the fraction 
of blue galaxies increases with projected distance, out to 
^ 1.5 Mpc* (as well as increasing dramatically with red- 
shift). Kodama & Bower (2000) confirm this strong trend 
in their recent photometric study of the intermediate red- 
shift clusters from the CNOC survey (Yee et al. 1996). Such 
trends of average colour with environment (i.e. radius or lo- 
cal density) could be due to (i) a systematic, intrinsic cluster 
colour gradient in one or both of the dominant elliptical and 
SO populations of the CMR with environment; or (ii) be pri- 
marily driven by the morphology-density relation (Dressier 
1980), where the late- type to early-type galaxies ratio is 
a strongly decreasing function of local galaxy density, or 
equivalently an increasing function of clustocentric distance 
(Whitmore & Gilmore 1991). Thus the observed trend may 
simply reflect evolution in the morphological mix, rather 
than representing changes in the colours of the galaxies ly- 
ing on the CMR with environment (as we discuss below, 
variations in the morphological mix within the CMR may 
produce more subtle changes in colour with environment). 

To reliably differentiate between these options we need 
to compare the properties of morphologically-selected (or, 
failing, that colour-selected) early-type galaxies in the out- 
skirts of clusters with similar galaxies in the high density 
cores. In this way we can determine if the blueing trend oc- 
curs in the passive cluster population or simply reflects an 
increasing proportion of late-type galaxies. 

Unfortunately, very few of the existing studies of clus- 
ters have the field of view, photometric precision (the ex- 



* Throughout this work values of Hq = 50 km s ^ Mpc ^ and 
qo = 0.5 have been adopted. 



pected colour shifts are ;c 0.10 mags) or colour information 
necessary to investigate the possible variation with environ- 
ment of the colour of the galaxies on the CMR. One of the 
few relevant studies is the recent investigation by Terlevich 
ct al. (2001) of the (U — V) colours of ~ 100 morphologically- 
classified early-type galaxies across a 1.5° field within the 
Coma cluster [z = 0.023). They identify a trend for the 
colours of early-type galaxies to be systematically bluer at a 
fixed luminosity outside the core of the cluster. They inter- 
pret this as a result of age differences and suggest that there 
is a 2 Gyr difference between the mean luminosity-weighted 
ages of the stellar populations of galaxies in the core and 
those outside 0.5 Mpc. 

At higher redshifts, deep multi-colour imaging with 
Hubble Space Telescope {HST) of distant clusters has suf- 
ficient sensitivity to uncover the signs of the variation in 
the CMR with environment, as well as providing detailed 
morphological information which can bo used to remove the 
effects of changes in the morphological mix from the analysis 
(Ellis et al. 1997). However, the limited field of view of the 
WFPC2 camera means that these HST-hased surveys are 
restricted to the inner regions of clusters. For example, van 
Dokkum et al. (1998) discuss the radial dependence of the 
CMR for the distant cluster MS 1358-1-62 aX z = 0.33 within 
a ~ 1.6 Mpc region. They find no radial dependence of the 
zero-point of the CMR, or the scatter around it, for bulge- 
dominated galaxies at any radius. They do find, however, an 
increasing scatter around the CMR for early-type disk galax- 
ies in the outer parts of their survey region (van Dokkum 
et al. 2000). They propose that this trend is caused by in- 
falling galaxies at the outskirts of the cluster, but the ab- 
sence of detailed morphological classifications for the early- 
type disk galaxies still leaves open the possibility that their 
result may again be due to subtle changes resulting from the 
morphology-density relation. 

Spectroscopy provides a much more sensitive tool to in- 
vestigate the variation in the ages of the stellar populations 
of galaxies with environment, as well as providing unambigu- 
ous membership information for galaxies at large projected 
distances from the cluster centre. Spectral signatures of past 
activity can remain measurable for ~2 Gyrs after the star 
formation has ceased, by which time the photometric signa- 
tures may have become undetectable (e.g. Terlevich et al. 
1999). Spectroscopic surveys of both local and distant clus- 
ters have uncovered precisely these signatures: identifying 
a class of galaxies with spectral features (enhanced Balmer 
absorption lines) indicative of recently ended star formation 
(Couch & Sharpies 1987; Caldwell et al. 1993). The galax- 
ies are classified as a-fk or k-l-a based on the strength of 
the Balmer absorption lines (Dressier et al. 1999), equiv- 
alent terms are E-l-A, PSG and HDS (Couch & Sharpies 
1987). Global gradients in these spectral indices have been 
identified in galaxies from a sample of X-ray selected distant 
clusters (Balogh et al. 1999). However, the samples available 
from most spectroscopic studies are relatively modest and 
for the most part lack detailed morphological information. 
Thus there is a similar ambiguity to that in the analysis 
of the photometric gradient, due to the morphology-density 
relation, and it is difficult to interpret the origin of the ra- 
dial gradients observed in key spectroscopic emission line 
indicators (Poggianti et al. 1999). 

One survey which combines the photometric and spec- 
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Table 1. Details of the clusters from the LARCS sample used in this work. For each cluster we give the coordinates of the cluster 

centre, the redshift (z) and its X-ray luminosity, Lx, in the 0.1-2.4keV passband (Ebeling ct al. 1996). We also list the parameters of 
the biwcight fit to the CMR within 2Mpc of the centre of each of the clusters as plotted in Figure 1. The slope, colour at the fiducial 
magnitude and associated errors are taken from the median of 100 realizations of the background subtraction. These errors may be 
underestimated due to the way in which the CMR is fitted (see text for details). The apparent i?-band magnitude corresponding to our 
fiducial absolute magnitude of My = —21.8 is also given for each cluster. The colour of each CMR at a fiducial magnitude equivalent to 
Mv = —21.8 and the slope of the relationship are examined as a function of redshift in Figure 2. 



Cluster 




R.A. Dec. 

(J2000) 




z 


Lx 

(lO'i'icrgs-l) 


Slope 


(B — R)My = -1\.i, 


flMv = -21.8 


A 22 


00 


20 38.64 


-25 43 


19 


0.131 


5.31 


-0.054 ± 0.008 


1.85 + 0.08 


17.12 


A 550 


05 


52 51.84 


-21 03 


54 


0.125 


7.06 


-0.057 + 0.006 


2.05 + 0.07 


17.01 


A 1079 t 


10 


43 24.90 


-07 22 


45 


0.132 


<0.45 


-0.045 + 0.003 


1.80 + 0.05 


17.14 


A 1084 


10 


44 30.72 


-07 05 


02 


0.134 


7.42 


-0.051 + 0.009 


1.83 + 0.06 


17.17 


A 1285 


11 


30 20.64 


-14 34 


30 


0.106 


5.47 


-0.056 + 0.014 


1.87 + 0.13 


16.61 


A 1437 


12 


GO 25.44 


+03 21 


04 


0.133 


7.72 


-0.072 ± 0.005 


1.94 + 0.12 


17.15 


A 1650 


12 


58 41.76 


-01 45 


22 


0.084 


7.81 


-0.039 ± 0.003 


1.69 + 0.06 


16.06 


A 1651 


12 


59 24.00 


-04 11 


20 


0.084 


8.25 


-0.047 + 0.003 


1.63 + 0.08 


16.06 


A 1664 


13 


03 44.16 


-24 15 


22 


0.127 


5.36 


-0.060 + 0.015 


1.73 + 0.13 


17.05 


A 2055 


15 


18 41.28 


+06 12 


40 


0.102 


4.78 


-0.052 + 0.013 


1.70 + 0.08 


16.53 


A 3888 


22 


34 32.88 


-37 43 


59 


0.151 


14.52 


-0.097 + 0.017 


2.12 + 0.16 


17.47 



(f) Although not originally fulfilling the LARCS selection criteria, A 1079 resides in the field of A 1084 at a similar redshift to the 
latter and is thus included in this study. 



troscopic techniques is the study of A 2390 {z = 0.23) by 
Abraham et al. (1996). They combine a large spectroscopic 
survey covering a wide field (to provide spectral line indices 
and membership) and precise photometry. They demon- 
strate that the colours of the reddest spectroscopically- 
confirmed cluster members on the colour-magnitude relation 
become progressively bluer as a function of clustocentric dis- 
tance out to 5 Mpc, a trend which was also seen in Balmer 
line indices. They interpret these radial gradients as arising 
from a decrease in the mean luminosity- weighted age of the 
stellar populations in galaxies at larger projected distances 
from the cluster's centre. They suggest that this results from 
the truncation of star formation as galaxies and groups are 
accreted from the field onto the cluster outskirts, building up 
the cluster in an "onion-ring" fashion reminiscent of the hi- 
erarchical picture of cluster formation (Lacey & Cole, 1993; 
Cole et al. 2000). 

The impressive studies of Coma by Terlevich et al. 
(2001) and A 2390 by Abraham et al. (1996) illustrate the 
power of photometric comparisons of galaxies across a range 
of environments to address the questions posed earlier. Both 
studies, however, only cover a single cluster and so it is not 
clear how general their conclusions are, especially given the 
cluster-to-cluster differences which could arise from different 
evolutionary histories. 

This paper describes an analysis of the Las Campanas 
Observatory and Anglo-Australian Telescope Rich Cluster 
Survey (LARCS) to address this problem. LARCS is a long- 
term project to study a statistically-reliable sample of 21 of 
the most luminous X-ray clusters at intermediate redshifts 
(z = 0.07-0.16) in the southern hemisphere. We are mapping 
the photometric, spectroscopic and dynamical properties of 
galaxies in rich cluster environments at 2: ~ 0.1, tracing the 
variation in these properties from the high-density cluster 
cores out into the surrounding low-density field beyond the 
turn-around radius. For the most massive clusters at « ~ 0.1, 



the turn-around radius corresponds to roughly 1 degree or a 
10 Mpc radius (O'Hely et al. 1998) and therefore we have 
obtained panoramic CCD imaging covering 2-degree diam- 
eter fields, as well as spectroscopic coverage of these fields 
(Pimbblct et al. 2001; O'Hely 2000; O'Hely et al. 1998). The 
imaging comes from B and J?-band mosaics taken with the 
1-m Swope telescope at Las Campanas Observatory, while 
the spectroscopy comes from the subsequent follow-up with 
the 400-fibre 2dF multi-object spectrograph on the 3.9-m 
Anglo- Australian Telescope (AAT). 



The advantage of this sample for the present analysis is 
the homogeneity of the clusters and observations. This al- 
lows us to combine many clusters together to provide a large 
sample of cluster galaxies across a wide range in environ- 
ment. In particular, we are able to use statistical corrections 
to remove field contamination, rather than requiring spec- 
troscopic membership. We can thus investigate the colour 
of galaxies across a wide range in environments in clusters, 
spanning three orders of magnitude in projected galaxy den- 
sity from the core out to close to the turn-around radius at 
~ 8 Mpc. 



In §2, we describe the LARCS cluster sample and the 
precision of our photometric catalogues. In §3 we outline the 
statistical method used to correct for field contamination 
and the construction of colour magnitude diagrams. In §4, 
we present our results on the characteristics of the galaxy 
populations in the clusters. By creating a composite cluster 
we explore the variation in the CMR with radius and local 
galaxy density. We discuss these results in §5 and present 
our conclusions in §6. 
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Figure 1. The colour-magnitude diagrams for tlie 11 clusters used in this work. All galaxies statistically-identified as cluster members 
within a 2 Mpc radius of the cluster centre from one realization of the field subtraction are plotted as heavy points. The smaller points 
are the galaxies which were rejected as field by our statistical subtraction method. The biweighted fit to the OMR is shown as a solid 
line. The Icr uncertainty in the colour is shown by the flanking parallel dashed lines. The vertical line denotes the absolute magnitude 
limit of My = —20 used to deflne the galaxy sample used for the fit. 



2 CLUSTER SELECTION AND 
OBSERVATIONS 

The sample used in our analysis is from the LARCS survey 
(Pimbblet et al. 2001). LARCS is a survey of the proper- 
ties of an X-ray luminosity-limited sample of 21 clusters at 
2 = 0.07-0.16. The clusters in the sample are selected from 
the X-ray brightest Abell Clusters (XBACS) catalogue of 
Ebehng et al. (1996). XBACS comprises 242 Abell clusters 
detected in the ROSAT All-Sky Survey and is effectively a 
complete, volume-limited sample of the most X-ray lumi- 



nous clusters, Lx ^ 5 x 10**ergs~^, at « ~ 0.1. The 21 
clusters in the LARCS sample are a random subsample of 
the southern XBACS clusters with Lx > 3.7 x lO^'^ergs"^ 
at z = 0.07-0.16. 

Ten clusters from LARCS are analysed in this paper 
and we list them in Table 1 (an additional cluster, A 1079, is 
present in our observations of A 1084 and we have included 
it in our analysis). Those ten clusters possess the highest 
quality and most complete photometric observations from 
the full survey. The subset covers a representative range of 
the complete LARCS sample, including clusters which ap- 
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pear relaxed in their X-ray and galaxy distributions (e.g. 
A 1650), as well as more dynamically young, disturbed look- 
ing systems (e.g. A 1664 - which appears to have undergone 
a recent cluster-cluster merger; Edge, priv. comm.). 

High quality broad-band B- and _R-band CCD images 
of each cluster have been obtained at Las Campanas Ob- 
servatory using the 1-m Swope Telescope. More details of 
the observations, reduction and analysis of these data arc 
given by Pimbblet et al. (2001); here we briefly summarise 
the pertinent points. 

A 2°-diamctcr field around each cluster is imaged in 21 
over-lapping pointings to produce a mosaic of a region out to 
~ lOMpc radius at the survey's median redshift of 2: ~ 0.12. 
The images are reduced using standard tasks within IRAP 
and are then catalogued using the SExtractor package of 
Bertin & Arnouts (1995). Wc adopt "mag_best" from SEx- 
tractor as the estimate of the total magnitudes. To deter- 
mine colours for the sources we measure aperture photome- 
try within 4" diameter apertures (typically ~ 10 kpc at the 
cluster redshifts) on seeing-matched tiles using phot within 
IRAF. Photometric zeropoints arc computed using frequent 
observations of standard stars (typically over 100 standard 
stars per broad-band filter per night) from Landolt (1992) 
interspersed throughout the science observations. On pho- 
tometric nights, the variation in the resultant colour and 
extinction terms are all within la of each other and the 
final photometric accuracy is better than 0.03 mags. 

As some of the observations within a mosaic are un- 
dertaken in non-photomctric conditions it is necessary to 
calibrate these independently. This is achieved through the 
method developed by Glazebrook et al. (1994). We use in- 
formation on the photometry of objects in the overlapping 
regions of the tiles to calculate the offsets between pho- 
tometric and non-photometric mosaic tiles. Once corrected 
by Glazebrook's algorithm, the resultant median tile-to-tile 
variation in zero-points across the mosaic is < 0.006 mags 
(see Figure 1 of Pimbblet et al. 2001), as estimated from 
the scatter in the magnitude from the duplicated sources in 
the overlapping regions. The maximum deviations between 
tiles is ~ 0.015 mags. We expect that the photometric ze- 
ropoint errors of 0.03 mags dominate the systematic error 
budget and the tile-to-tile variation in zero-points are small 
by comparison. 

Finally, we apply corrections for galactic reddening 
based on Sclilcgcl ct al. (1998). The final internal magni- 
tude errors across the full mosaics are typically 0.03 mag 
and always less than 0.06 mags (Pimbblet et al. 2001). Such 
accuracy is essential to detect subtle photometric radial gra- 
dients in the colour-magnitude relation. The catalogues are 
typically 80 percent complete at a depth of i? ~ 22.0 and 
B ~ 23.0. 

Star/galaxy separation for the catalogues uses the ro- 
bust criteria described in Pimbblet ct al. (2001): FWHM > 
2.0" and CLASS_STAR< 0.1 from the SExtractor star-galaxy 
classifier. The resulting stellar contamination based on these 
criteria is estimated at < 3 percent (see Pimbblet et al. 
2001). To further check for misclassification of compact 
early-type cluster galaxies as stars we have constructed 
colour- magnitude diagrams for the star samples and con- 
firm that they exhibit no hint of a CMR at the relevant 
colour for the cluster members. 



3 CONSTRUCTION OF THE 

COLOUR-MAGNITUDE DIAGRAMS 

In this section we describe the construction of colour- 
magnitude diagrams from the LARCS catalogues, the use 
of the biweight method to fit the CMR and the spatial dis- 
tribution of the cluster galaxies (i.e. cluster morphology). 

3.1 Statistical field correction 

To examine the galaxy population of the clusters in the ab- 
sence of spectroscopy it is necessary to correct for the "field" 
contamination in a given cluster catalogue. The statistical 
subtraction method utilized in this study is described in de- 
tail in Appendix A. 

Briefiy, the field population is determined from the 
outer regions of the LARCS cluster mosaics at a radius well 
beyond 6Mpc (at a redshift of z = 0.12). The field regions 
from each cluster arc all examined by eye prior to being 
stacked to give the average field distribution. After examin- 
ing these data, we reject part of the field around Abell 1084 
due to the presence of Abell 1079. The field sample is not 
significantly biased to a lower moan density because of its re- 
moval. Once generated, the final stacked field population^ is 
then appropriately area-scaled to that of the cluster sample 
we wish to correct. 

The colour- magnitude diagrams for galaxies from the 
cluster sample and the final, scaled, stacked field sample 
are then compared. By direct comparison of corresponding 
regions on the colour-magnitude diagrams we assign each 
galaxy in the cluster sample a probability of being a field 
galaxy. Then, using a Monte Carlo method, we subtract off 
the field population based upon these probabilities. As de- 
scribed in Appendix A, however, problems can arise if the 
calculated probability exceeds 1.0. The solution utilized in 
this work is to expand the interval in colour and magnitude 
used to calculate the probability until it lies in the range 
0.0 < P(Ficki) < 1.0 (see Appendix A for a full discussion). 

In this work, the statistical Monte Carlo back- 
ground subtractions are realized 100 times and the colour- 
magnitude fits presented represent the median of those back- 
ground subtractions. 

3.2 Fitting the CMR 

We run the background subtraction algorithm (see Ap- 
pendix A) on each cluster to produce 100 realizations of 
the cluster's colour-magnitude diagram. In fitting the CMR 
we only use those galaxies whose absolute rest frame magni- 
tude is brighter than My = —20. This magnitude is chosen 
to correspond to the equivalent limit in Butcher & Oemler 
(1984). 

Each of the CMRs arc fitted using a robust biweight 
method (Beers et al. 1990) identical to that employed by 
Terlevich (1998). Briefiy, the biweight method is a statis- 
tically robust estimator that weights outlying points (such 
as the blue population) low in order to find the best fit 

t The final effective area of our field sample is ~ 13.1 square 
degrees. Each cluster mosaic covers ~ 3.1 square degrees on the 
sky. 
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line to the CMR. In line with Beers et al. (1990), we adopt 
the median absolute deviation as the best estimator of the 
biweight location and scale and a value of the weighting pa- 
rameter of w = 6, this results in zero-weight being given to 
galaxies more than away from the median relation with 
the weighting increasing quadratically closer to the median 
value (e.g. 95 per cent weighting at la). Since we seek to fit 
a model to these data, we must minimize the residuals of the 
data to the model. Therefore the set of points Xi = Yi—f{Xi) 
arc defined, where {Xi,Yi) are the data points and f{X) is 
the model (a straight line of the form f{X) = mX + c). The 
scale and location of the then iteratively minimized 

on the (rn, c) plane using the multidimensional downhill sim- 
plex method of Press et al. (1992). We apply this routine to 
the background-subtracted cluster sample from each of the 
100 Monte Carlo simulations. We adopt the median of the 
fit parameters from the 100 simulations as our best estimate 
of the CMR and use the scatter between the sinnilations to 
give the uncertainty in this fit. We note that our adoption 
of the errors from the Monte Carlo field correction may un- 
derestimate the uncertainty in the fits when the sample is 
dominated by the cluster population. Most of our analysis, 
however, lies in the regime where the field contribution is im- 
portant and so we retain these errors as the best estimates 
of the internal uncertainty in our fits. 

In Figure 1 we show the distribution of galaxies on 
colour-magnitude diagrams for one such realization of each 
of the clusters within a radial extent of r-p < 2 Mpc. De- 
tails of the median values of the fits to the 100 realisations 
for each cluster and their associated errors are presented in 
Table 1. 

All the clusters show colour-magnitude relations for the 
redder, early-type galaxy members. We see strong CMRs 
(more than 60 galaxies lying on the CMR within 2 Mpc 
of the cluster centre) in A 550, A 1285 and A 3888; whilst 
A 1079, A 1084 and A 1664 exhibit much weaker relations 
(less than 25 galaxies lie on the CMR within 2 Mpc of the 
cluster centre). All the CMR's exhibit a negative slope as 
expected from the colour-luminosity relation of early-type 
galaxies (Table 1, Figure 2). 

The colour at a fiducial magnitude of Mv = —21.8 
(equivalent to an L* galaxy) is also calculated from the re- 
sult of the biweight fit; {B — R)mv=-2i.8- These colours are 
presented in Figure 2 and tabulated along with the values 
of Mv = -21.8 in Table 1. The {B - R)mv=-2i.8 colours of 
galaxies on the CMRs in the clusters form a tight relation, 
median deviation of 0.08 mags, which steadily reddens with 
redshift due to the K-correction. We also show the expected 
variation with redshift in the observed colour of a galaxy 
which formed all its stars at very high redshift, z ^ 2. 
The reddening in colour at the fiducial magnitude we see 
is broadly in line with that expected for a population which 
formed the bulk of its stars at high redshift (Metcalfe et al. 
1991). The small scatter in the redshift trend implies that 
our cluster to cluster photometry is internally consistent. 

The statistical field correction we apply is designed to 
detect overdensities above the level of the "average" field 
on the colour-magnitude plane. Therefore, if there are other 
ovor-dcnsc regions in the fields of our clusters wo expect 
to see some secondary CMRs in the colour-magnitude di- 
agrams. Such CMRs are seen in the colour-magnitude diar 
grams of A 550, A 1285 and A 1651. As expected from the 
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Figure 2. a) Median observed {B — R) colour at My = —21.8 
calculated from 100 realizations of the background subtraction 
plotted versus cluster redshift (with no K-correction applied) . The 
solid line is the best fit to the data, with the parallel dashed 
lines being ±1(t. The dotted line illustrates the expected evolution 
in the galaxy colours of galaxies which formed all their stars at 
high redshift, z 2> 2. The clusters form a tight relation evolving 
redward with increasing redshift. b) The variation in slope 
of the CMR with redshift. There is a trend for steeper CMR's 
at higher redshifts, as expected from the effects of shifting K- 
corrections on the observed {B — R) colours of the galaxies along 
the CMR. We show a "No evolution" model which shows the 
behaviour expected just from the differential K-corrections. 



available volume, these CMRs are higher redshift clusters 
in the background of the observed clusters, for example, we 
identify the background cluster in A 1651 {z = 0.084) as 
A 1658 {z ^ 0.145). To reduce the influence of these systems 
on the biweight estimator used to fit the CMR of the LARCS 
cluster, we employ an appropriate red limit to ensure the fit 
converges on the desired CMR. 
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3.3 The spatial distribution of cluster galaxies 

We now use the CMRs to map the two dimensional galaxy 
distribution within the clusters. Assuming that the galax- 
ies on the CMR are associated with the cluster, we define 
a cluster galaxy to be within the la scatter of its fitted 
colour- magnitude relation (listed in Table 1). This criteria 
is applied to the whole sample of galaxies in each field. Hav- 
ing obtained their positions on the sky, we use a circular 
top-hat function with a 500 kpc smoothing length to gener- 
ate smoothed galaxy distribution maps for each cluster. The 
resultant distributions for the central 25' of each cluster are 
presented in Figure 3. 

The brightest galaxy in these X-ray luminous clusters 
is typically a giant elliptical and we identify these using 
the CMR and the morphology of their extended low surface 
brightness halo in our imaging data. For the more regular 
clusters, the position of this galaxy is generally in agreement 
with the peak of the distribution of colour-selected cluster 
members (Figure 3) and we therefore use the position of this 
galaxy to define the cluster centre. Our clusters demonstrate 
a broad range of morphological types ranging from very reg- 
ular, highly-concentrated examples (e.g. A 550), through to 
very diffuse systems (e.g. A 1664) whose central galaxy po- 
sition does not coincide with any of the main spatial over- 
densities seen in the colour-selected cluster members. 



4 ANALYSIS AND RESULTS 

In this section we describe the analysis of the properties of 
the cluster's CMRs as a function of radius and local galaxy 
density. 

4.1 Radial variation of the CMR 

We search for any radial variations in the CMR within the 
clusters by dividing the central 6Mpc diameter region of 
clusters into three anulii. We note that the virial radii for 
these X-ray luminous clusters should lie in the range 2-4 
Mpc (Carlberg et al. 1996). For each radial bin, a biwoightcd 
fit is calculated to the CMR and the colour at a fiducial 
magnitude (equivalent to My = —21.8) is obtained. We 
have also calculated the colour based on a fixed slope (using 
the fit to the CMR in the central 2 Mpc diameter annulus) 
for the field-subtracted cluster realizations and found that 
this method produces the same results as an unconstrained 
gradient fit, within 1 standard deviation (see Table 2). The 
result of those fits are presented in Table 2, whilst the colour- 
magnitude relations for the three radial bins in each cluster 
can be seen in Figure 4. 

In most of our clusters we observe a modest blueward 
shift in colour at the fiducial magnitude with radius by 
A{B - R) -0.09 ± 0.11 out to 3 Mpc (Table 2). We see 
no significant change in this trend with cluster redshift. Dis- 
crepant results arise in three clusters. In A 1084, the change 
in (B — R) colour from the central radial bin to the third is 
much more than A{B — 7?) ~ —0.10. We attribute this to 
the poor fit in the third radial bin (see Figure 4) . The third 
radial bin of A 3888 is similarly poorly fitted. In A 1664, no 
radial trend is apparent. 

However, we note that there are two shortcomings of 



this analysis. Firstly, the formal errors on the biweighted 
fits for individual clusters are large enough that the colour 
gradient trend in the individual CMRs is insignificant. Sec- 
ondly, the biases in the statistical correction of the colour- 
magnitude diagram make it increasingly unreliable at large 
radii (Appendix A). We must therefore adopt a different 
approach. 

To trace the variation in the colour of the CMR out to 
larger radii and better quantify its strength, we can exploit 
the homogeneity of our cluster sample and observations by 
combining all of our clusters together into a single composite 
system. This will improve the statistics for cluster members 
at large radii and provide a more robust measurement of the 
trends in the CMR. In doing so, however, we acknowledge 
that there may be aperture bias present in our photometry: 
we use a fixed photometric aperture and therefore cluster 
galaxies will be sampled to different physical radii as a func- 
tion of redshift. This effect may alter our determined CMR 
slopes and colour changes if there are strong colour gradi- 
ents present within our cluster galaxies. We estimate that 
at least 85 per cent of the total light falls into our large 4" 
apertures (~ 10 kpc at 2 ~ 0.12; see §2) in all but the largest 
of our galaxies. As the nature of this work is a differential 
radial analysis of the composite CMR, we consider such an 
effect to be small. 

To produce a composite cluster we have to transform 
the magnitudes, colours and positions of the galaxies in the 
individual clusters onto a common scale. The transforma- 
tion of the magnitudes uses the relative magnitudes of non- 
evolving early- type L* galaxies from Table 1, an assumption 
which is supported by the agreement between the no evo- 
lution predictions for colour evolution and the observations 
across the limited redshift range spanned by our sample. 
The apparent magnitudes of the galaxies can therefore be 
transformed to a median redshift of z = 0.12. To trans- 
form the colours of the galaxies on the CMR onto a uni- 
form scale we simply reduce the observed CMR's to provide 
colours relative to the fitted relation. These can then be 
easily transformed to z = 0.12. The slope we adopt for this 
transformation is assumed to be constant with radius, and 
equivalent to the measured values presented in Table 1. We 
emphasize that our slopes are determined in a totally ho- 
mogeneous manner for a homogeneous sample: two absolute 
requirements given the uncertainties. The most meaningful 
method to compare the radial positions of galaxies in the 
various clusters is to use their radius normalised to the virial 
radius of the clusters, i?vir. However, the weak dependence 
of the virial radius on the X-ray luminosity, Rvir oc L^^® 
(Babul et al. 2001), combined with a narrow X-ray luminos- 
ity range for our sample, means that i?vir is only expected 
to vary by 20 percent across the whole sample. Hence, for 
simplicity, the clusters are all scaled to a fixed metric size. 

The final step in our analysis of the composite cluster 
is to construct colour histograms down to the fiducial mag- 
nitude limit, Mv = — 20. These histograms are statistically 
corrected using the binned colour distribution from the field. 
We can use this approach as we are interested in the typical 
colour of the CMR, rather than its slope, so we no longer 
need to retain the magnitude information on the galaxies. 
The advantage of using colour histograms is that it circum- 
vents some of the concerns about the field correction tech- 
niques used on the colour-magnitude plane (Appendix A). 
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Figure 3. The two dimensional smoothed spatial distributions of the cluster galaxies lying within the scatter of the biweight fit to the 
CMR. A circular top-hat function with a 500 kpc smoothing length is used on the positions of galaxies to generate these maps. The 
lowest contour represents a surface density of cluster galaxies of 6.0Mpc~^, with each subsequent inward contour denoting an increase 
of 2.0Mpc~^. The crosses indicate the adopted centre of the clusters based on the position of the apparent brightest cluster member 
from the CMR. 
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Figure 4. (a) Colour-magnitude diagrams for A 22, A 550, A 1079, A 1084, A 1285 and A 1437 in 1-Mpc wide radial bins. Each diagram 
shows one of the 100 realizations made with the background subtraction algorithm. The solid line, where present, is the biweighted fit 
to the OMR; if absent, a biweighted fit could not be successfully made. The short dashed lines denote the combined errors associated 
with the Tp < 1 Mpc fit of each cluster. The longer dashed, vertical lines denote the My = —20 fiducial magnitude limit; only galaxies 
brighter than this arc considered in the biweighted fit. For most of the clusters, there is a trend for the OMR to evolve bluewards by 
A{B -R) 0.09 ±0.11 out to 3 Mpc. 



We create the colour histogram of the field sample in 
an identical manner to the cluster samiilc. An appropriately 
area-scaled region of the field is subtracted off the compos- 
ite cluster histogram in each radial bin. The result of this 
analysis is illustrated in Figure 5. 

Figure 5 shows that the red cluster members in the 
CMR are still visible as a peak in the colour histogram as far 
out as ~ 6 Mpc from the cluster core. Beyond this, the red 
galaxies become swamped by an increasingly blue cluster 
population and the noise from the statistical field correc- 



tion. To identify the peak in the CMR and quantify how 

it changes with radius we fit the colour histograms with a 
gaussiaii using a minimization method. To make the fit 
we restrict the colour range to {B — R) = 1.6-2.2. An initial 
estimate of the modal colour, its amplitude and the width 
of the distribution, ctpeak, is made and fed in as the in- 
put to the gaussian fitting routine. The routine outputs new 
estimates of these values which are iteratively returned to 
the fitting routine until convergence is achieved. The resul- 
tant values are checked by visually inspecting the fit. The 
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Figure 4. (b) As for Figure 4a, but with A 1650, A 1651, A 1664, A 2055 and A 3888. 



gaussian fits generate the colour of tlic CMR at tfie fiducial 
magnitude and a measure of the width of the CMR, together 
with associated errors, in each radial interval. The fits are 
presented in Table 3 and illustrated in Figure 5. The peak 
of the CMR in the combined cluster evolves bluewards with 
radius at a rate oi d{B - R) / drp = -0.022 ±0.004 (where rp 
is the projected radius from the cluster core), equivalent to 
A{B — R) ~ —0.11 ± 0.05 from the central radial bin out to 
6 Mpc. Beyond 6 Mpc the constraints on the colour gradient 
trend of the CMR become less reliable and so we limit our 
discussion to the region within 6 Mpc. 

4.2 The colour— local density relation 

As the clusters in our survey exhibit a wide range in their 
morphologies (see Figure ,3) , we seek a more general method 
of combining the galaxy samples from the different clusters, 
rather than simple radial-averaging, to investigate the en- 
vironmental differences within the composite sample. Local 
galaxy density provides a comparable measure of environ- 
ment across the sample and so we employ this instead. 



We calculate the local projected density for each galaxy 
in our clusters and use it to assign the galaxy to its relevant 
environmental bin. The local galaxy density, E, is estimated 
by finding the surface area on the sky occupied by that 
galaxy and its ten closest neighbours down to the fiducial 
magnitude limit of My = — 20. Due to field galaxy contami- 
nation, this local density will be over-estimated. We correct 
for this by subtracting off a constant density estimated from 
the median local density of our combined field sample. We 
divide our sample up into eight logarithmically-spaced bins 
covering three orders of magnitude in local galax;y density 
(Table 3). 

For each local density bin a colour histogram is created. 
As before, the clusters are transformed to a median red- 
shift of 2 = 0.12 and the slope of the CMR is corrected for. 
Since these colour histograms are not corrected for field con- 
tamination, it is necessary to subtract off an appropriately 
area-scaled colour distribution for our field sample (see the 
lower-right panel in Figure 5). Finding the area from which 
a restricted range of densities are drawn is non-trivial; a re- 
gion of high density will have a much smaller field correction 
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Table 2. Parameters of the biweight fit to the CMRs of the clusters, split into three equal radial bins as plotted in Figure 4. The slope 
and colours are those found using the biweight estimator fit to each radial bin. The values reported in the final column (fixed slope) 
uses the slope of the central bin in the particular cluster. The errors are the la variation from the 100 Monte Carlo simulations of the 
background subtraction. 



Cluster 


Radius 
(Mpc) 


Slope 


{B — R)mv = -21.8 


(-B — -R)mv=-21.8 
Fixed Slope 


A 22 


0-1 


-0.054 ± 0.008 


1.85 ± 0.07 


1.85 ± 0.07 




1 o 

i— z 


— U.liU it U.UiU 


1 Qv -1- n 1 o 
l.o 1 ± v.Lz 


1 7Q _L n OI 

i. / o it U.zi 




2-3 


-0.052 lb 0.007 


1.75 ±0.17 


1.73 lb 0.22 


A 550 


0-1 


-0.057 lb 0.001 


2.05 ± 0.07 


2.05 lb 0.07 




1 o 
L — Z 


— U.UOU it U.UUo 


o no _i_ n 1 r; 
Z.XJZ it U.lo 


o no -L n 1 yi 

Z.\JZ it U.-L4 




2-3 


-0.061 ±0.014 


1.94 ±0.35 


1.92 ± 0.38 


A 1079 


0-1 


-0.045 ± 0.003 


1.80 ±0.05 


1.80 ±0.05 




1 o 
1— z 


— U.lUi it U.UU4 


1 OO J- A r\'7 
L.oz ± U.U ( 


1 Td J- A nn 

i. / D it u.uy 




z— O 


— U.UOi it u.uuy 


1 Qo -U n on 

L.oZ it U.ZU 


1 Qc; _|_ n 
i.c50 it U.oo 


A 1084 


0-1 


-0.050 lb 0.009 


1.83 lb 0.06 


1.83 lb 0.06 




L—Z 


—u.uoz ± U.U14 


1 '7Q -i- r\ AO 
L. (o it U.4z 


1 TA _L A Q1 

i. / y it U.oi 




z—o 


n nc;'7 _i_ n no'7 


1 _L. n Kyi 

i.OO It U.04 


1 fin _L A 
i.DD it U.OD 


A 1285 


0-1 


-0.055 lb 0.007 


1.87 lb 0.12 


1.87 lb 0.14 




L—Z 

Z—o 


— U.UdU it U.Uio 


L.oO ± U.oo 


i.i$D it U.o4 


A 1437 


0-1 


-0.077 lb 0.005 


1.94 lb 0.10 


1.94 lb 0.12 




L—Z 
z—o 


n nRK _i_ n nnQ 
— U.UDo it U.UUo 


1 Q'7 _L n oi 
L.o 1 ± KJ.ZL 


1 Q'7 _L A OI 
L.O 1 ± KJ.ZL 


A 1650 


0-1 


-0.038 lb 0.001 


1.70 lb 0.06 


1.70 lb 0.06 




L—Z 


— U.UoU ± U.UUo 


1 R c: _i_ n no 
I. DO it U.Uo 


1 Rt^ _l_ A AQ 
l.DO it U.Uo 




Z—O 
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— U.U41 it U.Uli 


1 fi'i -U n n? 

i.Oo it U.U ( 


1 Rc; -I- A nfi 
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A 1651 


0-1 


-0.047 ± 0.002 


1.63 ± 0.07 


1.63 ± 0.08 




1-2 


-0.060 ± 0.011 


1.63 ±0.16 


1.64 ± 0.24 




2-3 


-0.048 ± 0.008 


1.59 ±0.15 


1.60 ±0.14 


A 1664 


0-1 


-0.062 ±0.012 


1.73 ±0.07 


1.73 ±0.07 




1-2 


-0.058 ±0.039 


1.76 ± 0.07 


1.76 ± 0.08 




2-3 


-0.076 ± 0.024 


1.76 ± 0.08 


1.76 ± 0.08 


A 2055 


0-1 


-0.054 ±0.007 


1.71 ±0.08 


1.71 ±0.08 




1-2 


-0.050 ±0.021 


1.67 ±0.15 


1.62 ±0.18 




2-3 


-0.054 ± 0.024 


1.67 ±0.26 


1.60 ± 0.28 


A 3888 


0-1 


-0.096 ± 0.015 


2.13 ±0.13 


2.12 ±0.13 




1-2 


-0.103 ±0.043 


2.08 ± 0.22 


2.04 ± 0.25 




2-3 


-0.193 ±0.013 


1.88 ± 0.35 


1.90 ± 0.36 



than a region of low density. We generate an adaptive map 
of the field region by binning up the positions of the field 
galaxies and taking the median density of each bin. The 
map is thresholded using density cuts equal to the eight 
logaritiimically-spaced bins in local galaxy density to gener- 
ate the area occupied by galaxies of that local galaxy density. 
This measure of the area is used to scale the colour distri- 
bution of the field sample and subtract it from each local 
density bin's colour histogram. 

The results of this analysis are presented in Figure 6. 
Our combined cluster evolves bluewards with decreasing lo- 
cal density at a rate of d(B-i?)/dlogio(E) = -0.076±0.009, 
equivalent to A{B - R) -0.20 ± 0.06 from the liigliest lo- 
cal density regime covered by our sample (logj^(,(E) > 2.5) 
to three orders of magnitude lower. The results from this 
analysis provide a comparable estimate of the shift in the 
CMR colour to that identified in §4.1 (for the typical clus- 



ter in our sample a local galaxy density of logiQ(S) ~ 0.5 
corresponds to a radius of ~ 2 Mpc). 



4.3 Width of the CMR 

From the histograms presented in Figures 5 and 6, we 
now examine how the width of the CMR varies with en- 
vironment. The gaussians fitted to these data in the figures 
are used to obtain crpeak which we use as an estimate of 
the CMR's width. These values are plotted in the upper 
panels of Figure 7, with tabulated values presented in Ta- 
ble 3. We find that the CMR peak does appear to broaden 
with increasing radius and decreasing local galaxy density 
by A(Tpeak(B — 7?) ~ 0.15 across the ranges studied, al- 
though there is considerable uncertainty in the individual 
measurements. 

To further investigate how the colour distribution of 
galaxies in the CMR varies with environment we also cal- 
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Figure 5. Colour histograms of the composite cluster from the center to a maximum radial extent of 8 Mpc. The field colour distribution 
is also pictured in the lower-right for reference (arbitrarily scaled). The best fit gaussians are shown for the CMR peaks in each radial 
bin. The solid vertical line is the colour of the CMR peak in the central bin. The dashed vertical line is the peak of each individual radial 
bin. A blueward shift in the peak of the CMR can be seen clearly even by the second radial interval, although the position of the CMR 
peak becomes more uncertain at large radii. The CMR evolves blueward at a rate of d(B — R)/drp = —0.022 ± 0.004 from the centre to 
the outer regions. Also note that the CMR appears to the broaden at larger radius. 



culate the colours of the 30"^ and 70*'^ percentiles from the 
red end of the colour distributions (Figures 5 & 6). The 
30*'' percentile represents the reddest members of the CMR 
peak, the typical colour of this population appears to remain 
nearly constant in colour (Figure 7) out to ~ 6 Mpc (i.e. the 
limit of the visibility of the CMR in Figure 5). Beyond 6 
Mpc, there is a reddening of this colour which we believe 
is due to CMRs of other higher redshift clusters. A simi- 
larly selected sample exhibits the same near constant trend 
in colour with local galaxy density. In contrast, the trend of 



the 70' percentile, representing the bluer members of the 
CMR, shows a strong shift to the blue with both radius and 
local galaxy density. 

Thus whilst the colour of the reddest members of the 
CMR appears to be relatively constant in different environ- 
ments, the width of the CMR peak broadens at larger radii 
and lower local galaxy density, reflecting an increasing tail of 
blue galaxies on the outskirts of the clusters (see Figures 5, 
6 & 7). 
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Figure 6. Colour histograms of the composite cluster colour magnitude relations in logarithmic density intervals after correcting for the 
slope in the CMR and field contamination. The best-fit gaussians are shown for the CMR peaks in each density bin. The solid vertical 
line shows the colour of the peak of the CMR in the highest density bin. The dashed vertical line is the peak in the particular density 
interval. The CMR peak evolves steadily blueward at a rate of d{B — iJ)/iilog^Q(S) = —0.076 ± 0.009 across three orders of magnitude 
in local galaxy density (S). 



4.4 Blue cluster galaxies 

We now compare the relative proportions of the red and blue 
galaxies in the clusters. To do this we quantify the fraction 
of non-CMR galaxies by deriving the Butcher-Oemler blue 
fraction (Butcher & Oemler, 1984), fs- The blue fraction is 
calculated for all of our clusters using the method described 
in Butcher & Oemler's (1984) landmark study. Briefly, to 
calculate /s we use only those galaxies whose magnitude is 
brighter than Mv = — 20 and lie within a radius of Rso'- the 
radius of a circle containing 30 percent of the cluster's pro- 



jected galaxy distribution. After background subtraction, 
the net number of galaxies in these populations are a mea- 
sure of cluster richness, denoted N^o- The blue fraction is 
the fraction of this population whose rest frame colour is 
A{B -V) ^ -0.2 bluer than the fitted CMR. The concen- 
tration of the cluster is defined as logjQ(i?6o/-R2o)'^ and is 



t The radius Ra is the radius of a circle containing a per cent of 
the cluster's projected galaxy distribution. 
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Figure 7. The .30*'' and TO"' percentiles of the colour histograms presented in Figures 5 & 6, measured relative to the red end of the 
distribution. These arc plotted in the two lower panels as a function of radius and local galaxy density (solid lines) together with the 
peak CMR colour (dashed line) taken from Figure 5. The 30"^ percentile of the distributions is near constant over at least the inner 
6 Mpc of our combined cluster whilst the peak shifts to the blue by —0.11 ± 0.05 mag. We suggest that the sharp redward trend seen at 
large radii may result from the presence of background clusters and groups which dominate the red wing at the largest radii. The upper 
panel displays the width of the gaussian, (Tpeak, fit to the colour histograms against radius. The width steadily increases by ~ 0.15 mag 
across the radial range. Similar behaviour is seen in the distributions when sorted in terms of local galaxy density. 



a measure of the cluster's central concentration (Butcher & 
Oemler 1978, 1984). 

We present in Table 4 the median value of fs for the 
clusters and the associated values of R-m), A^^o and the con- 
centration. Also noted are the morphologies of the clusters 
from Figure 3. 

We stress that these values of fs are likely to be slightly 
low due to the biases in the statistical background correction 
technique employed. Sec Appendix A for discussion on this 
point. There may also be biases arising from the field sample 
we use. Since they are taken at the very edge of the cluster 
they may be fractionally more dense than the 'average' field 
thus further underestimating /s . We consider this effect to 
be small and the statistical background correction to be the 
dominant source of any bias. 

There is considerable scatter in the value of /b in the 
clusters in our sample. In the majority of our clusters the 



blue fraction is negligible and, combined with the small 
range in X-ray luminosity and redshift covered by our sam- 
ple, it is difficult to distinguish any trends in these param- 
eters. There arc two clusters, however, which do have sig- 
nificant blue populations; A 1437 and A 3888. Both of these 
clusters show evidence for recent merger events (Edge, priv. 
comm.) and appear highly concentrated in X-ray. The large 
scatter in /s for the clusters in our homogeneous sample is 
surprising and suggests that this parameter is sensitive to 
short term events which affect the galaxy mix within the 
clusters. 



5 DISCUSSION 

We have examined variation of the CMR with environment; 
both with radius and with local galaxy density. We find that 
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Table 3. Peak colour and full width of the CMR's variation 
with radius and local galaxy density as estimated from the fitting 
gaussians to the colour distributions (Figures 5 & 6). 



Sample Peak (B - R)my=-21.s o-peak 



Radius (Mpc) 

0- 1 1.92 ±0.01 0.13 ±0.01 

1- 2 1.88 ±0.02 0.15 ±0.02 

2- 3 1.84 ±0.03 0.18 ±0.03 

3- 4 1.83 ±0.04 0.20 ±0.03 

4- 5 1.81 ± 0.05 0.23 ± 0.04 

5- 6 1.81 ± 0.05 0.19 ± 0.04 

6- 7 1.83 ±0.05 0.23 ±0.05 

7- 8 1.84 ±0.06 0.25 ±0.06 

Log Density (Mpc"^) 

>2.5 1.96 ±0.01 0.15 ±0.02 

2.0-2.5 1.88 ±0.02 0.15 ±0.02 

1.5-2.0 1.87 ±0.03 0.17 ±0.03 

1.0-1.5 1.85 ±0.04 0.15 ±0.04 

0.5-1.0 1.76 ±0.04 0.19 ±0.04 

0.0-0.5 1.77 ±0.05 0.22 ±0.05 

-0.5-0.0 1.79 ± 0.05 0.22 ± 0.06 

-1.0- -0.5 1.76 ±0.06 0.23 ±0.08 

< -1.0 1.78 ±0.07 0.23 ±0.11 



Table 4. Median values for the Butcher-Oemler blue fraction, Jb , 
taken from 100 realizations of the statistical background subtrac- 
tion. The errors on /g are one standard deviation from the me- 
dian. Values of R30, ^30 and the concentration, Cone, (Butcher 
& Oemler 1979) derived from the calculation of fg are also tab- 
ulated. The cluster's two-dimensional distribution of CMR mem- 
bers is also qualitatively described: Regular, R; Irregular, I. 



Cluster 


/b 


R30 (Mpc) 


N30 


Cone 


Morph. 


A 22 


0.02±0.02 


1.06 


44 


0.29 


R 


A 550 


0.04±0.01 


1.35 


96 


0.32 


R 


A 1079 
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0.29 
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0.28 
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A 2055 


0.06±0.03 


1.26 


36 


0.29 
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A 3888 


0.21±0.02 
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63 


0.36 
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the CMR peak evolves bluewards and broadens with both 
radius and local galaxy density (Figures 5, 6 & 7). Yet the 
colour of the reddest members of the CMR remains constant 
(Figure 7). Wc emphasize at the outset of our discussion 
that the effects of environment upon galaxy evolution (e.g. 
star formation suppression, triggering and the morphologi- 
cal evolution of galaxies) are still polemical issues in modern 
astronomy (Dressier et al. 1997; Lubin et al. 1998; Andreon 
1998; van Dokkum et al. 1998). In this section, we exam- 
ine possible causes for these effects with particular reference 
to the SO evolution proposed by the MORPHS team (e.g. 
Dressier et al. 1997; Poggianti et al. 1999). 



5.1 The effects of morphology on the CMR 

We now estimate the trends which should exist in our sam- 
ple due to the existance of a morphology-radius relation (T- 
R, Whitmore & Gilmore 1991; Whitmore et al. 1993) or a 
morphology-density relation (T-E, Dressier 1980). We adopt 
a null hypothesis that the colour gradient we observe in the 
CMR with environment is due to an increasing fraction of 
early-type spiral galaxies in the cluster sample, while the 
colours of the spheroidal galaxies (ellipticals and SOs) are un- 
affected by their environment. This hypothesis is suggested 
by the fact that we see a broadening of the colour- magnitude 
relation with radius, but a fixed red envelope (Figure 7). 

Although we really do not have detailed morphological 
classifications for our sample, we can use the distribution of 
galaxies in terms of concentration index (CI) and peak sur- 
face brightness (/Umax) to look for changes in the morpho- 
logical mix in different environments within the clusters. On 
the CI-pMAx plane, early-type galaxies populate the high- 
concentration and high surface brightness region, with later- 
types typically having lower concentrations and fainter peak 
surface brightnesses (Abraham et al. 1994; Pimbblet et al. 
2001). Comparing the distribution in CI-^max of galaxies 
lying on the CMR in the cluster core (vp < 1 Mpc) and out- 
skirts {vp > 4-6 Mpc) we also find evidence for a shift in the 
typical morphology - with the outer regions having a laxger 
proportion of lower-concentration, lower-surface brightness 
galaxies than the core. To quantify this difference, we use a 
two dimensional Kolmogorov-Smirnov (K-S) test (Fasano & 
Franceschini 1987) to show that the CI-/iMAx distribution 
for the rp < 1 Mpc and rp > 4-6 Mpc samples are unlikely to 
be drawn from the same parent population at 97.5 per cent 
confidence. Unfortunately, it is non-trivial to transform the 
observed shift in the CI-/LtMAx distribution between these 
environments into expected colour differences. 

There is another approach we can use to test our null 
hypothesis: we can evaluate the typical {B — R) colour of 
the galaxy distribution in the clusters based on Whitmore 
& Gilmore's (1991) radial morphological mix and Dressler's 
(1980) local-density morphological mix. We adopt the mor- 
phological mix in local rich clusters on the assumption that 
it is unlikely to change dramatically from « = to a = 0.12 
(Fasano et al. 2001; Dressier et al. 1997). Our limiting 
magnitude of My = —20 is also broadly comparable with 
Dressler's (1980) limit of My ~ -20.4. We use the {B - V) 
colours of the different morphological types from the Third 
Reference Catalogue (RC3, de Vaucouleurs et al. 1991) as 
illustrated in Roberts & Haynes (1994) and assume a no- 
evolution K-correction based on the SEDs in King & El- 
lis (1985). We split the sample into broad morphological 
types corresponding to those employed by Dressier (1980) 
and Whitmore & Gilmore (1991): E, SO and S+Irr. The 
{B — V) colours of E and SO galaxies from the RC3 are 
{B — V) = 0.90 and 0.89 respectively. These correspond to 
observed {B - R) colours at 2 = 0.12 of {B - R) = 1.94 
for E's and {B - R) = 1.91 for SOs. The third morphologi- 
cal class is "S+Irr" : comprising spiral and irregular galaxies. 
The use of this class in our analysis, however, is problematic 
due to its broadness. The CMR will not contain many of the 
later-type spirals (e.g. Sc and Sd) as their mean colours lie 
significantly bluewards of the CMR defined by elliptical and 
SO galaxies. 
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Table 5. Relative changes in the colours of the cluster popu- 
lation a.t z = 0.12, compared to the galaxy population in the 
cluster core, assuming the morphological mixes from Whitmore 
& Gilmore (1991) and Dressier (1980). A blueward shift of -0.07 
mags in {B — R) is predicted from the centre out to 4 Mpc and 
—0.08 mags over a two orders of magnitude change in local density 
from the core. 



Sample E SO Sa/Sab A{B - R)My = -2i.s 



Radius (Mpc) 

0- 1 0.27 0.67 0.06 0.00 

1- 2 0.25 0.60 0.15 -0.04 

2- 3 0.23 0.59 0.18 -0.05 

3- 4 0.21 0.56 0.23 -0.07 

Log Density (Mpc~^) 

>2.5 0.45 0.55 0.00 0.00 

2.0-2.5 0.43 0.53 0.04 -0.02 

1.5-2.0 0.36 0.55 0.09 -0.04 

1.0-1.5 0.29 0.56 0.15 -0.07 

0.-5 1.0 0.21 0..59 0.20 -0.09 

0.0 O..") 0.20 0.57 n.2:i -0.10 



To examine which morphological typos from the 
"S+Irr" class are likely to inhabit the CMR wc estimate 
how blue a galaxy has to be before it falls outside of the 
CMR peak and hence does not affect our gaussian fitting to 
the colour distributions (Figures 5 & 6). Wc estimate that 
galaxies as blue as (B — R) ~ 1.5-1.6 will affect the fitting of 
the CMR in the outskirts of the clusters; this corresponds to 
{B — V) ^ 0.7 at z = or the colour of an Sab galaxy from 
RC3. At z = 0.12, an Sa galaxy will have a mean colour of 
{B-R) = 1.55 while Sab's have (B - R) = 1.49. Hence, wc 
expect that it is only the Sa/Sab class from the spiral pop- 
ulation which is likely to infiuence the colour of the CMR 
based on our fitting procedure. We therefore use Sa/Sab as 
our third morphological class, instead of the broader "S-fIrr" 
class. Before wc calculate the expected colour of the CMR 
peak we re-normalize the T-R and T-E morphological mix- 
tures to account for using only Sa/Sab galaxies. This again 
is done using the RC3, from which we estimate the frac- 
tion of Sa/Sab galaxies within the "S+Irr" class and then 
rc-scalo the relative proportions of the three morphological 
classes (Table 5). 

The gradient in mean {B — R) colour predicted from 
the variation in morphological mix seen in the T-R relation 
(Table 5) is d{B - R)/drp = -0.017. This is equivalent to a 
change of A{B — R) = —0.07 out to 4 Mpc which is roughly 
75 percent of that observed in the combined LARCS clusters 
{d(B-R)/drp = -0.022 ±0.004). The gradient in the CMR 
colour resulting from the variation in the morphological mix 
seen in the T-E relation is d(B- J?)/(ilogio(E) = 0.033. This 
is equivalent to a change from logjg(E) > 2.5 to logj^Q(E) = 
0.0 of A{B - R) = -0.10. which is around 50 percent of 
that seen over this density range in the composite of the 
LARCS clusters {d(B - i?)/dlogi(,(E) = -0.076 ± 0.009). 
We find similar predicted trends at 99 per cent confidence 
if we restrict ourselves to only the high Lx clusters from 
Dressier (1980). 



5.2 Interpretation 

Based on the results from the previous section wc conclude 
that the morphology-radius and morphology-density rela- 
tions seen in local clusters can account for only around half 
of the apparent change in galaxy colour with radius and 
local galaxy density seen in our sample. 

We estimate therefore that there is an intrinsic gradi- 
ent corrected for morphological differences in the colours 
of early-typo galaxies at a fixed luminosity of d(B — 
7?)/c?logj^Q(E) — —0.011 from the cores of typical rich clus- 
ters at « ~ 0.12 out to environments with galaxy densities of 
logiQ(E) ~ (equivalent to a radius of 2 Mpc). The intrinsic 
variation with clustocentric radius gives a somewhat weaker 
trend (Table 5): d(B - R)/drp = -0.005. 

Radial gradients in the colours of cluster galEixies have 
been previously reported by Abraham et al. (1996) and Ter- 
Icvich et al. (2001) for two rich clusters straddling the red- 
shift range studied fioro. Tficsc two measurements arc not 
exactly comparable as the Terlevich et al. analysis focuses 
on the morphologically classified E and SO galaxies, whereas 
Abraham et al. investigate a colour-selected sample around 
the CMR, similar to the approach used here. Nevertheless, 
as Terlevich et al. stress, the core of Coma is completely 
dominated by E and SO galaxies and so their results would 
not change significantly if they considered the whole popu- 
lation. 

Abraham et al. (1996) find radial gradient in {g — r) 
colour of d{g — r)/dlogio(rp) = —0.079 out to 5 Mpc in 
A 2390 {z = 0.23), this is not corrected for any variations 
due to a T — R relation. We estimate there is --^ 20 percent 
uncertainty in this gradient. Terlevich et al. (2001) report 
d{g - r)/dlogio(rp) = -0.024 ± 0.005 for the Coma clus- 
ter, transformed to the same restframe wavebands as the 
observations A 2390. To compare our result with these we 
similarly transform the {B — R) colours of galaxies in our 
composite z = 0.12 cluster to the wavelengths equivalent to 
observing a cluster in (gr — r) at 2 = 0.23. The result of this 
transformation is a radial gradient of d{g — r)/dlogi()(rp) = 
—0.061 ± 0.011 (uncorrected for the T — R relation). Hence, 
the radial gradient of the colour of the CMR at z = 0.12 
(a look-back time of 2 Gyrs) is intermediate between those 
found by Abraham et al. (1996) for A 2390 and Terlevich 
et al. (2001) for the Coma cluster at look-back times of 
3.5 Gyr and 0.5 Gyrs respectively. Correcting the higher red- 
shift measurements for the effects of the morphology-radius 
relation still results in steeper gradients at higher redshifts. 

Both Abraham et al. (1996) and Terlevich et al. (2001) 
attribute the colour gradients they observe in the clusters 
to age trends (Terlevich et al. 2001 use detailed arguments 
to show that the variation in colour cannot be due to dust 
spread through the core of the cluster) . Using the models for 
a single burst stellar population at an age of 7 Gyr (equiva- 
lent to the time elapsed between a formation epoch at 2: > 2 
and 2 — 0.12) presented in Bower et al. (1998), we find 
d{B-R)/dt ~ 0.03 magGyr"^ (Kodama & Arimoto, 1997). 
This implies that the galaxies on the CMR in the outer re- 
gions of our clusters have luminosity-wcigiited stellar popu- 
lations which arc approximately 3 Gyrs younger than those 
in the cores (r^ < IMpc). This age difference is compara- 
ble to that found by Terlevich et al. (2001) who predict that 
early-type galaxies in the outskirts of Coma have stellar pop- 
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ulations which are 2 Gyr younger than those in galaxies 
in the central regions. 

We also propose that the colour gradients observed in 
our sample represent variations in the ages of the stellar 
populations in the cluster galaxies arising from diffences in 
their star formation histories (Abraham et al. 1996; Small et 
al. 1998, 2001; Kodama & Bower 2000). The bulk of galaxies 
on the CMR are expected to be SO's (Table 5; Dressier 1980; 
Fasano et al. 2001) and thus changes in the colours of this 
morphological class must play a major role in causing the 
observed colour gradient (Abraham et al. 1996). However, 
we also see evidence for a proportion of old, evolved galaxies 
whose colours are constant with environment and define the 
red wing of the CMR (Figure 7). We suggest that the ma- 
jority of these galaxies are probably ellipticals, which show 
modest environmental variation in their properties (Treu et 
al. 1999). Thus a mix of ellipticals with constant colours 
and a dominant population of SO galaxies, whose colours 
vary with environment, can reconcile the main results from 
our survey. 

Why should galaxies which have SO morphologies, or at 
least the colours of SO's at the present-day, vary with en- 
vironment? Although still controversial, there is increasing 
evidence that the SO population in rich clusters has come 
into being relatively recently. The proportion of SO galaxies 
in rich clusters appears to increase rapidly to the present- 
day (Dressier et al. 1997; Fasano et al. 2001), possibly as 
the result of the transformation of star-forming disk galax- 
ies which are accreted by the clusters from the surrounding 
field (Poggianti et al. 1999; Kodama & Small 2001). The fact 
that we can see differences in the properties of cluster galax- 
ies out at least as far as 6 Mpc from the cluster core suggests 
that the process which may suppress the star formation of 
infalling galaxies can operate in low density environments 
(Fasano et al. 2000; Bekki et al. 2001; Carlberg et al. 2001; 
Couch et al. 2001; Kodama et al. 2001). One mechanism for 
causing this transformation is the suppression of star forma- 
tion due to the removal of the gas reservoirs as galaxies enter 
the cluster (Poggianti et al. 1999; Balogh, Navarro & Morris 

2000) . This mechanism may lead directly to the formation 
of an SO, or other processes may be involved (Balogh et al. 

2001) . 

If this scenario is correct then the accretion and trans- 
formation of spiral galaxies from the surrounding field will 
lead to a gradual build-up of the SO population within the 
cluster. Some morphological studies of the early-type galax- 
ies in cluster cores at 2 ~ 0.5 imply that the majority of the 
SO population in these regions today are transformed in the 
6 Gyrs between z — 0.5 and the present-day (e.g. Poggianti 
et al. 1999). On the outskirts of the clusters we expect a 
higher proportion of recently arrived cluster members, re- 
fiecting the accretion history of the cluster (Abraham et al. 
1996). Thus the SO galaxies in these regions will show the 
strongest signatures of their past star formation activity, 
exhibiting the youngest luminosity-weighted ages for their 
stellar populations and hence the bluest colours. These two 
effects will combine to produce a radial colour gradient in 
the same sense as that seen in our data. After accounting 
for the varying morphological mix with environment, it ap- 
pears that galaxies on the outskirts of these clusters have 
luminosity-weighted ages which axe ~ 3 Gyrs younger than 
those in the core. It would be interesting to quantitatively 



compare this observation with theoretical predictions to test 
models for the infall history and transformation of galaxies 
in high density environments (Diaferio et al. 2001). 



5.3 Caveats 

Our analysis and interpretation is based upon a number 
of assumptions that require consideration and should be 
improved upon in future work. Perhaps the largest source 
of uncertainty is the background subtraction method em- 
ployed. It is apparent from Figure 1 that there are a num- 
ber of galaxies, which are rejected by the background sub- 
traction algorithm, but which lie within the CMR (e.g. 
Abell 1650). Conversely, there must also be galaxies which 
are defined as cluster members, but which are in all likeli- 
hood interlopers. We re-emphasize that the background cor- 
rection is purely a statistical method. There axe two clear 
ways in which this method could be improved. Firstly, we 
could use a larger sample of dedicated field observations to 
improve the statistics. Secondly, we could use spectroscopic 
observations of the clusters to provide a more robust method 
to determine cluster membership. 

During our analysis we combine our clusters into 
a single composite cluster to examine the radial depen- 
dance of the CMR. We note that the richer clusters (i.e. 
Abell 550, 1285 and 3888) have a greater weight in the com- 
posite cluster than the poorer ones (i.e. Abell 1079, 1084 
and 1664). Although cluster-to-cluster variations are present 
within our sample, we emphasize that LARCS represents a 
homogeneously X-ray selected sample of clusters and there- 
fore such cluster-to-cluster variations should be minimized. 
Using more clusters to create the combined cluster would be 
beneficial. 

In our discussion we compare our data to that of 
Dressier (1980) to remove the T — T, relation from our 
clusters. We note that the LARCS dataset does not con- 
tain sufficiently high resolution imaging required morpho- 
logically classify galaxies. Hence we are unable to derive 
a T — S relation for our clusters and we rely upon the 
morphology-density analysis performed by Dressier (1980) 
for other clusters. In using these data, we note that the 
Dressier (1980) sample is constructed with an inhomoge- 
neous local {z < 0.06) sample of clusters and may, therefore, 
not be directly comparable to our study (although limiting 
Dressler's sample only to high Lx clusters gives the same re- 
sults). To remedy this situation would require sub-arcsecond 
panoramic imaging of our clusters to provide the necessary 
morphological information. 



6 CONCLUSIONS 

We have analysed eleven clusters with precise photometric 
observations from the LARCS survey to trace the variation 
in the colours of evolved galaxies into the environmental 
transition region between clusters and the field. 

• All of the clusters show colour-magnitude relations; 
many of these are strong and can be traced well beyond 
the central 1 Mpc, out to unprecedented clustocentric radii 
coverage of ~ 8 Mpc. 



© 0000 RAS, MNRAS 000, 000-000 



18 Pimbblet et al. 



• We find that the characteristic colour and slope of the 
CMR of galaxies in the cluster cores change with cluster 
redshift across the range z = 0.07-0.16 covered by our sam- 
ple. The variations in the OMRs are compatible with those 
expected from galaxy populations which formed the bulk of 
their stars at high redshifts, z ^ 2. 

• The majority of clusters in the survey exhibit a vari- 
ation with radius in the colour of the CMR at a fixed lu- 
minosity. When combined into a single composite cluster, 
we observe a radial colour gradient of d{B — R)/drp = 
-0.022 ± 0.004 from the cluster centre to ~ 6Mpc. 

• The composite cluster formed from our sample also 
shows a colour-local-donsity relation for the galaxies on the 
CMR. We identify a colour gradient in the modal colours of 
galaxies on the CMR amounting to d{B — 7?)/dlogjQ(S) = 
—0.076 ± 0.009 across three orders of magnitude in local 
galaxy density from the high-density core of the cluster to 
the outskirts. 

• We calculate the expected change in the colour of the 
CMR with environment due to the changing morphological 

mix. We estimate that approximately 50 to 70 per cent of 
the radial and density colour gradient is due to an increase in 
contamination by quiescent, early-type spiral galaxies falling 
close to the CMR. We suggest that the CMR consists of 
a passively evolving population (elliptical and evolved SO 
galaxies) plus recent additions (morphologically young SO 
and recently star- forming spiral galaxies). 

• Our analysis of eleven clusters supports the previous 
claims of radial colour variations in individual rich clusters 
at z = and z = 0.23. We have also shown that similar 
trends arc visible when using local galaxy density as a more 
general mesisure of galaxy environment. Taken together with 
previous studies, we suggest that these environmental trends 
in the colours of galaxies at a fixed luminosity most likely re- 
flect differences in the luminosity- weighted ages of the galax- 
ies in different environments. If interpreted purely as a dif- 
ference in ages and accounting for differences in the mor- 
phological mix, then the gradient observed in our sample 
suggests that the luminosity-weighted ages of the dominant 
galaxy population within the CMR at 6 Mpc from the clus- 
ter core are ~3 Gyrs younger than those residing in the core. 

• The results of our photometric analysis can provide 
a readily testable, quantitative prediction. If the observed 
colour gradient simply reflects age differences in passive 
stellar populations then we would expect that the typical 
H6a + H'yA line strength of galaxies lying on the CMR 
should increase by ~ 4A when moving from the cluster core 
to the outskirts of these systems (Terlevich et al. 1999). 

• Further studies are urgently needed to investigate the 
properties of galaxies in the lower density environments 
which connect rich clusters to the surrounding field, partic- 
ularly at higher redshifts. This region is essential for inves- 
tigating the impact of environment upon the star formation 
and morphologies of galaxies. 



This is the second paper in a scries based upon the 
LARCS survey. In the next paper we will examine 2dF spec- 
troscopy for several of the clusters from the survey (O'Hely 
et al. in prep.). 
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APPENDIX A: STATISTICAL CORRECTION. 

To examine the cluster galaxy population in the absence 
of spectroscopy, it is necessary to statistically subtract the 



field population from a given cluster plus field population. 
Typically, to obtain the cluster population one subtracts 
off an area and density-scaled portion of the defined field 
population, Fpop, from the defined cluster plus field popu- 
lation, C + Fpop. Both the -Fpop and the C + fpop arc then 
binned onto a sensibly gridded colour-magnitude diagram. 
A straight forward subtraction is then performed to calcu- 
late the probability of a galaxy in a given colour-magnitude 
grid position of being a field galaxy. Thus: 



P(Field)coi,mag = 



N(Field) 



col, mag 



X A 



N(Cluster -|- Field)coi,i 



(1) 



where A is an area- and density-dependant scaling factor. 
For each galaxy in C + Fpop , we proceed to generate a ran- 
dom rmmber between 0.0 and 1.0 and compare it to the 
galaxies' value of P(Field)coi,mag in order to determine if 
the given galaxy should be in the resultant cluster popu- 
lation. Problems in this method occur when directly sub- 
tracting the populations on the color-magnitude plane pro- 
duces a negative number of galaxies in the cluster popu- 
lation, Cpopn, for a particular grid position. A resultant 
P(Field)coi,mag can thus be greater than 1.0. Therefore to 
utilize a background subtraction method to analyze a clus- 
ter population in the absence of spectroscopy, it is necessary 
to solve this "negative galaxy" problem. A typical method 
is to normalize all of the probabilities to the largest value 
of P(Field)coi,mag. Clearly such a solution has drawbacks as 
the resultant P(Field)coi,mag distribution no longer reflects 
the appropriately scaled input data and has actually been 
scaled away from its true value in many locations on the 
color-magnitude plane for the sake of one (or a few) anoma- 
lously low counts in C -I- i^pop for a particular grid position. 

The method adopted in this study is similar to that 
used by Smail et al. (1998) and Kodama & Bower (2001, 
KBOl). If, when gridding up the colour-magnitude dia- 
grams, we produce a negative galaxy in the computation 
of P(Field)coi,mag, we solve it by increasing the grid size 
for that particular grid position. Thus the grid position 
(col, mag) is increased to include its nearest neighbours 
(col : col + 1, mag : mag + 1) and then P(Field)coi,mag is 
recalculated using this new 2x2 grid. If this is not enough 
to cure the negative galaxy problem, we increase the grid 
size for that particular grid position by one further step to 
(col — 1 : col -|- 1, mag — 1 : mag -|- 1) and so on until the 
field-corrected galaxy count exceeds zero. 

In Figure 8 we compare the resultant P(Field)coi,mag 
distribution using this method and KBOl. The grid position 
(3, 2) has a fluctuation in C -|- Fpop so that the number of 
galaxies in the corresponding grid position in Fpop makes the 
resultant P(Field)3,2 > 1-0. By normalizing to the anoma- 
lous value, all the rest of the probabilities are reduced from 
their original values and we are still left with a high proba- 
bility at the original grid position. The KBOl method sub- 
tracts the excess probability (in this case 0.2) and distributes 
it evenly between its neighbours, making the original grid 
position P(Field)3,2 = 1-0. Using our new method, we re- 
sample grid positions (3 : 4, 2 : 3) to make one new, larger 
grid position and calculate the new probability for this. 

It is found that the negative galaxy problem exists in 
all clusters examined in this study. The number of grid- 
expansions required to cure the problem are typically around 
ten 2x2 expansions and a small number of larger order 
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Figure 8. A comparison of the statistical correction mctliods on the colour-magnitude plane. The galaxies have been binned into 
a colour-magnitude grid and the calculated P(Field) is shown for each grid position, [top left]: An example distribution of P(Field) 
containing a grid position with P(Field) > 1.0 - the negative galaxy problem, [top right]: Normalizing all the points to the anomalous 
grid position retains the original anomalous point and artificially decreases the probabilities in the remainder of the grid, [bottom left]: 
The approach used by KBOl: the excess 0.2 probability is distributed equally into neighbouring grid cells, [bottom right]: Adaptive-mesh 
method used in this work - the lowest possible number of grid positions are used to recalculate the P(Field) distribution and no obvious 
fluctuation remains. 



expansions. The expansions are typically situated away from 
the CMR, preferentially where the C + Fpop is sparse. 

As with all background subtraction methods, this tech- 
nique has advantages and drawbacks. Whilst it preserves the 
original probability distributions better than other methods, 
it also has a tendency to "smear" the resultant distribution. 
To illustrate this one can imagine a situation where one grid 
position has a P(Field) > 1 whilst its neighbours have very 
small values of P(Field). Upon expansion of the grid, the 
probabilities of being a field galaxy significantly increase for 
those grid positions that initially had a very low P(Field). 
Therefore we effectively remove galaxies from adjoining grid 
positions to make up for the one grid position which had 



P(Ficld) > 1. It is because of this that wo may expect that 
statistics on sparsely populated areas of a colour-magnitude 
diagram to be adversely affected (e.g. it is likely that the 
blue population, and hence /b, is underestimated). Unlike 
previous methods, this new technique minimizes the number 
of grid positions affected in this manner. 
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